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We prove that a system of discrete 2D in-plane dipoles with four possible orientations, interacting 
via a 3D dipole-dipole interaction plus a nearest neighbor ferromagnetic term, has periodic striped 
ground states. As the strength of the ferromagnetic term is increased, the size of the stripes in 
the ground state increases, becoming infinite, i.e., giving a ferromagentic ground state, when the 
ferromagentic interaction exceeds a certain critical value. We also give a rigorous proof of the 
reorientation transition in the ground state of a 2D system of discrete dipoles with six possible 
orientations, interacting via a 3D dipole-dipole interaction plus a nearest neighbor antiferromagnetic 
term. As the strength of the antiferromagnetic term is increased the ground state flips from being 
striped and in-plane to being staggered and out-of-plane. An example of a rotator model with a 
sinusoidal ground state is also discussed. 



1. INTRODUCTION 

Recent advances in film growth techniques and in the experimental control of spin-spin inter- 
actions have revived interest in the low temperature physics of thin films [5, 8, 15, 25, 30, 31, 
35, 36, 38, 41, 42, 44]. These quasi-2D systems show a wide range of ordering effects includ- 
ing formation of striped states, reorientation transitions, bubble formation in strong magnetic 
fields, etc. [3, 26, 39]. The origins of these phenomena are, in many cases, traced to competition 
between short ranged exchange (ferromagnetic) interactions, favoring a homogeneous ordered 
state, and the long ranged dipole-dipole interaction, which opposes such ordering on the scale 
of the whole sample. The present theoretical understanding of these phenomena is based on a 
combination of variational methods and a variety of approximations, e.g., mean- field and spin- 
wave theory [14, 22, 26, 28, 29]. The comparison between the predictions of these approximate 
methods and the results of MonteCarlo simulations are often difficult, because of the slow relax- 
ation dynamics associated with the long-range nature of the dipole-dipole interactions [14, 44]. 
It would clearly be desirable to have more rigorous results about the spontaneous formation of 
such patterns. In a previous paper [20] we began to investigate these questions by means of a 
spin-block reflection-positivity method which, combined with apriori estimates on the Peierls' 
contours, allowed us to: 

(i) Describe the zero temperature phase diagram of a ID Ising model with nearest neighbor 
ferromagnetic and long range, reflection positive, antiferromagnetic interactions. These include 
power-law type interactions, such as dipolar-like interactions. We proved in particular the 
existence of a sequence of phase transitions between periodic states with longer and longer 
periods as the nearest neighbor ferromagnetic exchange strength J was increased. 

(ii) Derive upper and lower bounds on the ground state energy of a class of 2D Ising models 
with similar competing interactions, which agreed within exponentially small terms in J with 
the energy of the striped state. 

In this paper we prove for some models of dipole systems with discrete orientations, on a 
2D lattice, that their ground states display periodic striped order. As in the ID case the stripe 
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size increases with the strength of the nearest neighbor exchange interaction, becoming infinite, 
i.e., giving a ferromagentic ground state, when the ferromagentic interaction exceeds a certain 
critical value. The proof is again based on a combination of reflection-positivity and Pcierls' 
estimates and on an exact reduction of our 2D model to the ID Ising model studied in [20] . The 
analysis takes explicit account of the tensorial nature of the 3D dipole-dipole interaction that, 
in the absence of any short range exchange, tends to produce order in the form of polarized 
columns (or rows) of aligned spins, with alternating polarization. For dipoles oriented along 
four possible directions at each site the ground state shown in Fig.l is 4- fold degenerate (spin 
reversal and 90° rotations) and these are the only ground states. See [14] for a more detailed 
description of pure dipole states and [18] for some proofs. 



The order described in Fig.l is induced by the fact that two dipoles with the same orientation 
attract if their axes are parallel to their relative position vector, and they repel if both their axes 
are perpendicular to their relative position vector: in other words the dipole-dipole interaction, 
although overall antiferromagnetic (in the sense that it prefers total spin equal to zero) is, 
roughly speaking, ferromagnetic (FM) in one direction and antiferromagnetic in the orthogonal 
direction. It is then possible to show that, in the presence of an additional nearest neighbor 
ferromagnetic exchange, the FM order will persist in one direction. In the orthogonal direction 
the column-column interaction can then be effectively described in terms of a ID model, which 
can be treated by the methods of [20]. 

A similar method allows us to investigate the zero temperature phase diagram of the 2D 
discrete dipole model in the presence of a nearest neighbor mitiferroiriagncdic exchange, even in 
the case that the dipoles are allowed to orient along six different directions (four in-plane and two 
out-of-plane). In this case we prove the so-called reorientation transition [1, 14], consisting in a 
flip from an in- plane ground state, like the one depicted in Fig.l, to an out-of-plane staggered 
state, as the strength of the antiferromagnetic exchange interaction is increased. 

The effects of the interplay between the tensorial dipole-dipole interaction and a short- 
range exchange interaction is still far from being understood in general systems of continuous 
(Hciscnberg-like) spins, or in general anisotropic models with different assumptions on the al- 
lowed spin orientations and/or on the short range exchange interactions. It is worth remarking 
that in these classes of models the naive procedure of minimizing the Fourier transform of the 
pair interaction not only is wrong (because it neglects the local constraints coming from the re- 
quirement that the norm of the spin vector at each site is equal to 1) but generally leads to very 
bad estimates on the ground state energy and on the size of the characteristic zero temperature 
patterns. However it is interesting to note that, for an 0{n) spin model, n > 2, with a scalar 



t 
t 
t 
t 
t 
t 
t 
t 



I 
I 
I 
I 
I 
I 
I 
I 



t 
t 
t 
t 
t 
t 
t 
t 



t 
t 
t 
t 
t 
t 
t 
t 



FIG. 1: A ground state of the pure dipole-dipole system in 2D. 
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interaction whose Fourier transform admits a non trivial minimum at ko 7^ 0, it is relatively 
simple to show existence of "soft" striped, i.e., sinusoidal ID spin wave, ground states. This 
was shown via a very general argument by Nussinov [33] for classical 0{n) spins. This general 
result is independent of the dimension and of any reflection-positivity of the interactions. It is 
applied here to the concrete example of rotators in two dimensions interacting with a nearest 
neighbor ferromagnetic interaction and a weak long range scalar interaction, i.e., a scalar 
0{n = 2) model with competing long range interactions, imitating the decay properties of the 
real dipolc-dipolc potential. The fact that the solution to this scalar isotropic problem is easy 
adds further motivation to the study of the harder and more exciting case of real dipole systems. 
(Note that even the existence of the thermodynamic limit in an external field, which is expected 
to be shape dependent, is unproven for a 3D dipole system). 

The rest of this paper is organized as follows. In Section 2 we introduce the 2D spin model 
with discrete orientations and state our results about its zero temperature phase diagram. This 
includes existence of striped order (in the case of a short range ferromagnetic exchange) and 
existence of a reorientation transition from an in-plane to an out-of-plane ordered state (in the 
case of an antiferromagnetic exchange). In Section 3 we discuss an example of a 2D rotator 
model with long-range scalar interactions whose ground staters arc; given by sinusoidal ID spin 
waves ( "soft stripes" ) . In Sections 4 to 7 we present the details of the proof for the 2D spin 
model with discrete orientations, both for the case of ferromagnetic and of antiferromagnetic 
exchange. In Section 8 we summarize some aspects of the conjectured positive-temperature 
behavior of dipole systems with competing interactions. The rigorous analysis of the phase 
diagram at positive temperatures is presently beyond our reach. 

Some results and proofs in the present paper rely on those of the previous paper [20]. We 
have discovered a minor technical error in the proofs of Theorems 1 and 2 of [20] which was 
caused by overlooking several exponentially small terms of the form e~'^^ which came from the 
use of periodic boundary conditions on a ring of length L. These errors can be repaired, but we 
found that everything can be done more easily and clearly using open boundary conditions on 
the line rather than periodic boundary conditions on the circle. This improved methodology, 
which might be independently interesting for future work, is given here in Appendix A. 



2. 2D DISCRETE DIPOLES: THE MODEL AND THE MAIN RESULTS 

In this section we introduce the 2D model with discrete orientations and state the main 
results on the structure of its zero temperature phase diagram. We shall first discuss the stripe 
formation phenomenon in the presence of a nearest neighbor ferromagnetic interaction and 
then the reorientation transition phenomenon in the presence of an antiferromagnetic exchange 
interaction. 



A. The ferromagnetic case 

Let A C be a simple cubic 2D torus of side 2L and let S'x, x e A, be an in plane unit 
vector with components {5'x}i=i,2- We shall assume that can only be oriented along the two 
coordinate directions of A {i.e., that and can take values {—1, 0, 1}, with (S^)"^ + (S^)'^ = 
1). We shall denote by fl^ the corresponding spin configuration space and, for later convenience, 
we shall define ^1 = {S^ e : S'i = 0, Vx G A} and nj[ = {|a G ^^a : 5'x = 0, Vx g A} to be 
the subspaces of vertical and horizontal spin configurations. The Hamiltonian is of the form: 

2 

i,j = lx,yGA <x,y>GA 

where, denoting the Yukawa potential by l^(x) = e~^l''l |x]~^, the interaction matrix W(x) is 
of the dipole form given by: 

Wi, (x) = ^ {-d^dj)Y,{ic + 2nL) , X ^ (2.2) 
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and Wij{0) = ^n^o{—didj)Y^{2nL). The second sum in (2.1) runs over pairs of nearest 
neighbor sites in A and the constants J and A will be assumed nonnegative. The A term is 
inserted in H to discourage neighboring spins from having orthogonal polarizations. It has the 
effect of encouraging stripes, but this term alone cannot create stripes. Without the J term 
the ground state would be as in figure 1. 

One of the main result of this section concerns the zero temperature phase diagram of model 
(2.1) with e = 0, that is in the case that the long range dipole-dipole interaction is the "real" 
dipolar one. It is summarized in the following Theorem. 

Theorem 1. Let e = and < J < Jo(0), where 

There exists \o{J) such that, if X > Ao(J), then the specific ground state energy of (2.1) in the 
thermodynamic limit is given by: 

lim ^Eo{K)=mme{h) (2.4) 

|A|^oo |A| he^+ 

where e{h) lim|A|^oo \M~^-^per{-^) and Epel{A) is the energy of a periodic configuration with 
either vertical or horizontal stripes all of size h and alternate magnetization. If the side of A is 
divisible by the optimal period (i.e., by 2h*{ J), with h*{ J) the minimizer of the r.h.s. of (2.4)) 
then the only ground states are the periodic configuration with either vertical or horizontal 
stripes all of size h*(J) and alternate magnetization. 

If X = +00, i.e., if we consider model (2.2) restricted to fi](uf2^, then the same conclusions 
are valid for all J >0. 

Remarks. 1) For any J < Jo(0), the ground state has non trivial stripes of finite size, 
that is the minimizer h*{J) defined in the Theorem is finite, and h*{J) diverges logarithmically 
as J tends to Jo{0). In the hard core (A = +oo) case, Theorem 1 implies that for any J > Jo(0) 
we have ft* (J) = +oo and the ground state is ferromagnetic. 

2) The constant Ao(J) introduced in the Theorem is proportional to h*{J), i.e., it is of the 
form Ao(J) = Ch*{J), with C independent of J. In particular Ao(J) diverges at J = Jo(0). 
It is not clear whether this divergence is an artifact of our proof and whether one should 
expect the same result to be valid for smaller values of A (in particular for A = 0). For 
A = 0, we tried to look for states with all four possible orientations and energy smaller 
than the one given by (2.4), but we did not succeed. It could very well be that the same re- 
sult is actually valid for any A > 0, but unfortunately we don't know how to prove or disprove it. 

In the case that £ > 0, we can extend the previous result to all values of J > and finite A. 
The result is summarized in Theorem 2. 

Theorem 2. The conclusions of Theorem 1 are also valid if: e > 0, J > and X > 
maxjCe — J, 0}, where 

c.=co„s.{i>'?/i 1^^;/^: (2.5) 



Remarks. 1) Similarly to the e = case, the optimal stripe size h*{J) is finite if and 
only if J < Jo(e), where 

with 



a(m, k,e) = y/ in'^m'^ + P + , (2.7) 



5 



The constant in the r.h.s. of (2.5) is chosen such that is larger than Jo(£), f^or aU e > 0. 
Note that the conclusions of Theorems 1 and 2 are also valid under the assumptions; e > 0, 
J > and A > const/i*(J). 

2) For small £, the bound on the value of A above which we get the striped state is proportional 
to |loge| and in particular diverges at £ = 0. Also in this case, it is not clear whether this 
divergence is just an artifact of our proof and whether one should expect the same result to be 
valid for smaller values of A. 

In summary, we prove that the ground state of (2.1) is a periodic array of stripes if the 
A term is large enough, and more precisely: in the dipole case, £ = 0, if A is larger than a 
constant proportional to the period of the A = +oo ground state; in the Yukawa case, £ > 0, 
if A is larger than a constant depending on the Yukawa mass (and diverging logarithmically as 
it goes to zero). The proof of Theorems 1 and 2 in the A = +oo case, described in Section 4 
below, is based on an exact reduction to an effective ID model, to be treated by the methods of 
[20]: as a byproduct of the proof, we also get explicit bounds for the energies of excited states. 
The extension to finite A, both in the £ = and £ > case, is based on a Peierls' estimate 
on the energy of a droplet of horizontal spins surrounded by vertical spins. The proof of the 
Peierls' estimate is very simple under the assumptions of Theorem 2, that is in the presence of 
exponential decay with mass e and of a large penalty A (large non uniformly in the mass £); 
this proof, together with the necessary definitions of Peierls' contours and droplets, is described 
in Section 5. Extending the result to the case £ = and A finite {i.e., the case considered in 
Theorem 1) is highly non-trivial: in fact, since the long range potential decays as the third 
power of the distance, the naive dimensional estimate on the dipole energy of a droplet of size ^ 
decreases as — ^log£. So in order to exclude the presence of droplets of spins with the "wrong" 
orientation we have to use cancellations, that is we have to prove the presence of screening in 
the ground state. The proof is given in Section 6. 



B. The antiferromagnetic case 

Let us now consider the case that the <Sx have six possible orientations, four in-plane and 
two out-of-plane, i.e., S^, and Si can take values {-1, 0, 1}, with (5^)^ + {S^)"^ + (5^)2 = l. 
Their interaction is given by the Hamiltonian (2.1), with x still located at the sites of the 2D 
torus A considered above, with J < and £ = 0. Then, 

Theorem 3. There exists an absolute constant Aq > such that, if X > max{Ao — |-/|,0}, 
then the specific ground state energy is given by: 

lim i^Eo(A)=min{ei,eo-2|J|} (2.8) 

|A|^oo |A| 

with ei the specific energy of the planar antiferromagnetic state described in Fig.l and eo — 2| J| 
the specific energy of the out-of-plane staggered state. One has eo > e\. For \ J\ < (cq — ei)/2 
and A large enough the planar antiferrom,agnetic states (i.e., the one described in Fig.l plus 
those obtained from it by translation and/or by a 90° rotation) are the only ground states. For 
\J\ > (eo — ei)/2 and A large enough the out-of-plane Neel states are the only ground states. 
For \J\ = (eo — ei)/2 and A large enough the out-of-plane Neel state and the in-plane state 
described in Fig.l are the only ground states. 



Theorem 3 is proven in Sec. 7. Its proof goes along the lines of the proof of Theorem 2. 
In other words, we first consider the restricted problem with the dipoles all constrained to 
be parallel (either in- or out-of-plane), and determine the ground states using a reflection 
positivity argument. Then we prove that for A large enough the picture doesn't change, using a 
Peierls' argument. Note that now the critical value Ao is an absolute constant, independent of 
J: the Theorem can be extended to cover the £ > case, in which case Aq is also independent 
of £. This is due to the fact that the reference states we now need to consider (i.e., the 
ground states of the restricted problem) are always antiferromagnetic with period two, and 
this induces a strong screening which effectively corresponds to a faster decay of interactions, 
at least as the inverse distance to the power four. 
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3. 0{n) SPINS WITH SCALAR INTERACTIONS 

In this section we describe an example of a spin system with dipole-type interactions that 
displays periodic striped ground states. The construction of the ground state, based on an 
observation about 0{n) spin models with scalar interactions, was originally given by Nussinov 
[33]. This shows that striped states are, very naturally, the ground states for a large class 
of 0{n) models with n > 2 and competing interactions. For this class of models the local 
constraint on the norm of the spin at each site is not strong enough to invalidate the naive 
procedure for determining the ground states based on minimization of the Fourier transform of 
the interaction. 

Let us consider an 0{n) model, n > 2, on a simple cubic c?-dimensional torus A c Z*^ of side 
L with d> 1 and Hamiltonian: 

// = 5^J(x-y)5.-SV (3.1) 

where J(x) is a function of the distance |x| only and the 5x are n-dimensional unit vectors, 
n>2. The Hamiltonian (3.1) can be rewritten in Fourier space as follows: 

n 

^ = E E •^(l^)^k-^-k (3-2) 

i=l keXij:, 

where D^, = {k = 27rL-% , n e [0, n Z'^} and 



1= ^ Sy--"^ , J(k) = ^ J(x)e-'' (3.3) 



k 

Proposition 1 (Nussinov). // ko is a minimizer for j(k) and g G 0(n) is an n x n 
orthogonal matrix, then all states of the form 

5^9.ko = g . with ^ ( cos(kjj . x) , sin(ko • x) , , • • • , 0) , (3.4) 

are ground states of model (3.2)-(3.3). 

Proof. By (3.2) and the normalization condition ^ I^^P = |A| we see that the groimd 
state energy cannot be smaller than |A| mink -^(k); on the other hand the states (3.4) all have 
energy equal to |A| mink -/(k) and they satisiy the normalization l^xp = 1, so this proves that 
they are ground states. ■ 

Note that if ko ^ 0, the ground states {S^^'''°}xeA are "striped", in the sense that they are 
ID sinusoidal spin waves in the direction ko/|ko| with wavelength 27r/|ko|; we call these soft 
stripes. Under additional conditions on the geometry of the set of minima of J(k) one could 
actually prove that such states are the only possible ground states, but we will not investigate 
this question in the greatest possible generality. Instead, as an illustration, we will discuss an 
explicit example where these ideas can be used to infer that all ground states are soft stripes. 
A similar analysis can be performed to show that the previous remark also applies to the case 
of a 3D 0(n) model, n > 2, with spins interacting via a short range ferromagnetic interaction 
plus a positive long range Coulombic 1/r interaction, i.e., to the class of models considered in 
[11]. 

Example. Let kdl? he a simple cubic 2D torus of side L and let H he defined as 

<x,y> mez2 ^'.ySA I I 

where the first sum ranges over nearest neighbor sites of A and the 's are 2D unit vectors. 
Then there exists a constant c > such that, if e < cJ and L is large enough, all the ground 
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states of H are soft striped states. 

Remark. This result, as well as the general remark above, depends crucially on the 
fact that the interaction is scalar and isotropic, i.e., for each pair (x, y), the coupling depends 
only on S'x • 5'y. In the case of anisotropic interactions, the ground state may look very 
different, and will be in general very difficult to identify. As an illustrative example, take the 
case of soft scalar spins, i.e., scalar spins constrained to satisfy \S]^\ < 1, interacting via a 
pair potential J(x) as in (3.1), with the further property that J(0) = 0. This example can be 
viewed as an extreme anisotropic case, where all couplings between the i-th components of the 
spins, i ^ 1, have been switched off. In this case the ground state will be of the Ising type; this 
can be shown in the following way. The Hamiltonian is linear in S^., for each x e A, and it 
has to be minimized midcr the constraint that \S^\ < 1, for all x € A. By the aforementioned 
linearity, the minimum will be attained at the boundary, that is = ±1, for all x e A. The 
problem then reduces to determining the ground state of an Ising model with both long range 
antiferromagnetic and nearest neighbor ferromagnetic interaction, that is a very difficult and, 
in many respects, open problem, see [20, 29]. For intermediate values of the anisotropy the 
problem is even more difficult and it is unclear whether the transition from a sinusoidal spin 
wave state to an Ising-like state as the couplings between the i-th components of the spins, 
i 7^ 1, are decreased will be a sharp transition or rather a continuous one. 

Proof of the Example. The Hamiltonian (3.5), up to an overall constant, can be rewritten 
as in (3.2), with J(k) = 2JJ2i{^ - cosfe) + eEx^o e*''''|x|-3. Note that J^^^o e^^^lxl'^ can 
be conveniently rewritten as ^(k) — 1 + X^pg27rz2 / c?xe'(P+''^^(x^ + 1)~^/^, where ^(k) is the 
Fourier transform of (7 (x) = (1 — 5x) [|x|~^ — (x^ + l)^''/^] (note that g(x) goes to zero as |x|~^ 
as |x| oo). Using the fact that / dxe'(P+'<^)''(x2 + l)-3/2 = 27re-IP+''l (see [21]), we find that 

2 

J(k) =2J^(l-cosfci)+27re^e-|P+''l+£3(k)-£ (3.6) 

i=l pjto 

where g(k) is twice differentiable and even in k. An elementary study of the minima of 
(3.6) shows that for e/J small they are located within 0{e/.JY from the points (±7r£/J, 0), 
(0, ±7re/J). Moreover the minima have the following property (that we shall call non- 
degeneracy): there are no two distinct (unordered) pairs of minimizing vectors (p, q), (p',q')! 
such that p + q = p' + q'; this is because the minima all lie on a surface of strictly positive 
curvature. This implies that all ground states are "striped" , that is they are all of the form 
(3.4), with stripes that arc all (almost) horizontal or vertical. In fact one can show that no 
state obtained as a superposition of different minimizing modes can be a ground state: this is 
proven by using the normalization j^xj^ = 1 and the non-degeneracy condition. For a similar 
discussion, see [33]. ■ 

4. 2D DISCRETE DIPOLES: THE A = +oo CASE 

As mentioned in Section 2 we shall first prove Theorems 1, 2 and 3 in the case A = +oo, that 
is in the case that the spins are either all vertical or all horizontal or, possibly, all out-of-planc. 
Then we will show that, if A is large enough, the spins in the gromid state configurations will 
in fact be either all vertical or all horizontal. We will do this first for the ferromagnetic plus 
Yukawa case (Section 5), then for the ferromagnetic plus dipole case (Section 6) and finally for 
the antiferromagnetic plus dipole case (Section 7). 

A. The ferromagnetic case. 

Let us assume here that £, J > in (2.1) and (2.2) and the spins are oriented along four 
possible directions, as described in Section 2 A. Let us denote by iJoo the Hamiltonian (2.1) 
with the hard core interaction corresponding to A = +oo. By definition, Hoo is the restriction 
of (2.1) to U f2^, ie.to the space of configurations with spins either all vertical or all 
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horizontal. Without loss of generality, let us consider the case that all spins are vertical, that 
is 5x = (0,(Tx), = ±1, Vx G A. In the following we shall alternatively refer to the spins as 
"up or down" spins or "plus and minus" spins, with interchangeable meaning. The key remark 
is that in this case the system is reflection positive with respect to ferromagnetic reflections 
in horizontal lines. The relevant reflection symmetry is defined as follows. Let tt be a pair of 
horizontal lines midway between two lattice rows which bisect the torus A of side 2L into two 
pieces A_|_ and A_ of equal size. Let r denote reflection of sites with respect to tt. Clearly 
rA_ = A+. We define 

0a^ = ar^ (4.1) 

For any function i^({crx}xGA), we shall define the reflected function 9F as OF({(Tx}xeA) = 
[F({0(Tx}xeA)]*- Note that if depends only on the spins in A+, then F_ = will depend 
only on the spins in A_. We shall say that is reflection positive (RP) with respect to reflec- 
tions in horizontal planes if it can be written in the form Hao = C+ (x)6C+ (x)dp{x) 
for a positive measure dp{x), with _ff+,C+(a;) depending only on the spins in A+. In our case 
this representation can be achieved by defining as the interaction of the spins in A+ among 
themselves, H- = 9H+ as the interaction of the spins in A+ among themselves and by suitably 
rewriting the interaction term X^xeA+ SyeA_ o'xO'y(W22(x — y) — J5|x-y|,i), in the following 
way. 

Let X2 > and let us rewrite 
where in the last expression kj^ = (^1,^3). If Xq, > 2/2 and 82 = 



^ (TxO-yW22(x - y) = - ^ ^ tTxCTyi9|y'£(x - y + 2ni) = 

yeA_ yeA_ 

= - ^ ^ ( dk^ y'ki + £2 ^ (7x try e*'=i ) e - V^i (4.3) 

yeA_ 

that has clearly the correct structure — / C{x)9C{x)dp{x). Note that the exchange interaction 
between A+ and A_ admits a representation of the form — JX^i '^i^o'ij with Ui the spins in A+ 
at distance 1 from A_. This concludes the proof that is RP w.r.t. reflections in horizontal 
planes. 

Let TZj be the i tli row of A and, given a spin configuration ct^^ = {(TxjxeA, let Si{g_j^) be the 
spin configuration = restricted to TZi. By refiection positivity and, more specifically, by the 
chessboard estimate (see Theorem 4.1 in [16]), the energy of a;^ can be bounded below as 

2L 

Hoo{g_^) > ^ E ^oo({5i(a^), . . .,Si{a^)}) (4.4) 

1=1 

where {^^(a;^), . . . ,5i(CT^)} is the spin configuration in A obtained by repeating in every row 
the same spin configuration 5i(a^). Then the; problem of minimizing H^o among all possible 
spin configurations is reduced to the problem of minimizing among the spin configuration 
in which all spins in a given column have the same direction. Let us consider one such config- 
uration. If column xi has all spins pointing up, then we will label it by a^^ = +, if it has all 
spins pointing down we will label it by o-a;i = — . The energy of the configuration labelled by 
{(^x}x=i,...,2L is 



E = 2L 



2L 2L 
xjty=l x=l 



(4.5) 



with Wo = Y,n>i W22{ne2) and 

w^^{x -y) = E ^22 {{x-y + 2Lni)ei + 71262) 
nez^ 
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-E 



r dk2dks _ j^2 ^ik2n2^-\x-v+2Lm\y/kl+kl+e^ (-4 gx 

J ^/fc2 + fc2 + £2 2 V • ; 



where ei, 62 are the two coordinate unit vectors and in the second line we used the representation 
(4.2)-(4.3). Performing the summations over ni,n2 wc got 

^TT(^-y) = E /'^^^(^^ 1 _ ,-.L.im,k,) ' (4-7) 



mez 



where a(?n, fc,e) was defined in (2.7). This potential decays exponentially at large distances 
and is reflection positive with respect to the antiferromagnetic reflections. The ground states 
of this onc-dimcnsional model can be studied by the same methods of [20]. As mentioned in 
the Introduction, the proofs of Theorems 1 and 2 in [20] actually contains a minor technical 
error, which was caused by overlooking exponentially small terms which came from periodic 
boundary conditions. While fixing the error, we discovered that the same results can be proven 
by a simpler and more elegant method, which is presented here in the Appendix. The result is 
that the ground states for the one dimensional system described by the effective Hamiltonian 
(4.5) consist of periodic arrays of blocks of alternating sign, all of size h, with h the positive 
integer minimizing the function: 

2J 2 /-^ „ 47r2m2 e'^^"'*^'^) ha{m,k,e) 
e(/i =wi-2J+ — -- V / dk— — ^tanh ^ ' ' (4.8 



mez 



where wi = wo + X]n>i '^TT ('^)- Note that e{h) is the specific energy of a periodic configurations 
with blocks all of the same size and alternating sign. An elementary study of the behavior of 
e{h) shows that if J > Jo, with Jo defined as in (2.6), then the ground state is ferromagnetic. 
If < J < Jo then the ground state is non trivial and h is the integer part of the solution to 



Jo-J = 2> / dk 5^ 2 ^ o ' (4.9) 

^J-oo a(m,fc,£) • (1 - e-«('"''='^)) • (l + e-''«('"''=>^)) 



In terms of the original dipole system this shows that for J > Jo the ground state is ferromag- 
netic, while for < J < Jo the ground state is striped, with stripes either all horizontal or all 
vertical. The stripes have alternating orientation and their thickness varies from 1 to 00 as J 
is increased from to Jo- 

The results in [20] and in the Appendix also imply a bound on the energy of any given 2D 

state CT^ different from the striped ground state ct^. From now on, let us assume for simplicity 
that the minimization problem m.uif^^^+ e{h) is solved by a unique h* with the property that 
e{h*) = min^g2+ e{h) = min;,gR+ e{h). Note that this is not a generic property: in general 
min/j£2+ c(^-) 7^ min/j£j{+ e(/i) and for some special values of J it could even happen that the 
minimization problem on Z is solved by two consecutive values h*,h* + 1. However the 
general case can be treated in a way completely analogous to the one described below, at the 
price of slightly more cumbersome notation [45]. So let us assume that e(h*) = min^gg+ e{h) = 
min/j£,{+ e{h) and let cj — i min/i^/i» (e(/i) — e{h*)). Let us think a^^ as a collection of ID 
configurations corresponding to the configurations Si{a^) on different rows (here i = 1, . . . , 2L 
is the index labelling different rows). Moreover, let us think each Si{a/^) as a collection of blocks 
B of size hs and alternate magnetization. Eq.(23) in [20] implies that the energy of ct^ can be 
bounded from below as: 

H^{aj,)>\A\e{h*)+ ^ 2cj/ib (4.10) 

So for any block of wrong size ft-s we pay a penalty cjhs and for any block of optimal size 
h* we don't pay a priori any penalty. We can actually improve the estimate above in the case 
that a block of optimal size on a given row is next to a non-optimal block of size Hb on the 
same row: in this case we can use the proof of Eq.(26) in [20] to infer that for each such pair 
we pay a penalty 2dj{h* + hs), with dj = \ min/i^ft,* (e({/i, h*}) — e{h*)), and e{{h, h*}) is the 
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specific energy of a periodic configuration with blocks of sizes {. . . ,h,h, h*, h* , h, h, h*, h*, . . .). 
Combining (4.10) with the refinement we just discussed we find: 

H^{af,)>\k\e{h*)+ J2 ''jhB+ J2 dj{hB,+hB,) (4.11) 

B-.hBT^h" <Bi,B2>: 

where the second sum runs over pairs of nearest neighbor horizontal blocks, such that one of 
the two blocks in the pair is of optimal length. A straightforward computation, based on the 
explicit expression (4.8) of the specific energy, allows one to check that cj, dj > Ke~"^ , for 
suitable constants a, k > 0. This concludes the proof of Theorems 1 and 2 for the A = +oo 
case. 



B. The antiferromagnetic case. 



Let us now assume that £ > 0, J < and that the dipoles can have six possible orientations, 
as described in Section 2B. In the A = +oo case, the spins can only be oriented all parallel to 
each other, either out-of-plane or vertical in-plane or horizontal in-plane. 

If the spins are all oriented out-of-plane, then the Hamiltonian reduces to a long-range anti- 
ferromagnetic Ising Hamiltonian of the form: 

Hl= SlW^,{^-y)Sl + \J\ (4-12) 

x,yeA <x,y>eA 

with 5^ = ±1. As a consequence of the analysis in [17] (see Proof of Theorem 5.1) the ground 
state of the Hamiltonian (4.12) is the usual period-2 staggered state, for any value of | J| > 0. 
Its ground state energy is eo — 2| J|, with 



/ dk — — -— ^^TT > (4-13) 



where Wq^ = X^„>i(— l)"W^ii(«'e2) and a{m,k,e) was defined in (2.7). In order to prove 
Theorem 3 in the A = +oo case we need to compare this specific energy with the specific energy 
of the best possible in-plane spin configurations. 

So let us consider the case that the spins are all oriented in-plane; we can assume without loss 
of generality that they are all horizontal. (The choice of horizontal rather than vertical spins is 
made here in order to keep the definition of the pair tt of refiection planes the same as in the 
previous subsection, see the lines preceding (4.1).) In this case we define the variables CTx in 
such a way that /Sx = (cx, 0). Now the relevant reflection in horizontal planes, replacing (4.1), 
is Oa^ = —(Trx (note the minus sign!). By the chessboard estimate we reduce to an expression 
analogue to (4.5), given by 



E = 2L 



2L 2L 
xjty=l x=l 



(4.14) 



with ax now representing the direction (right or left) of the spin in (x,0), Wq^ the constant 
defined after (4.13) and 



g-\x-y\a{m,k,e) _|_ ^-(2L-\x-y\)a(m,k,e) 



w{x-y) = - Y / dka{m,k,s) _ 2La(m,k,e) • (4-15) 



For |J| = the groimd state of (4.14) is ferromagnetic (that is, in terms of the original 2D 
dipole model, the ground state is given by Fig.l - rotated by 90°). If | J| is increased then, 
using the same methods of [20], one can prove that the ground state has a sequence of transi- 
tions from the ferromagnetic state to periodic states of antiferromagnetic blocks of size h with 
alternating staggered polarization (of period h or 2h, depending whether h is odd or even). 
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As an illustration, the states with /i = 2,3,4 (with periods 4,3,8, respectively) are given by: 

( \ hH ), ( \ hH hH \ ) and ( 1 \ \ h + -H ). 

For any value of \ J\ the optimal period h*{\J\) is obtained as usual by minimizing over h the 
specific energy e"^^{h) of the states with AF blocks all of size h and alternating staggered 
magnetization. Here e^^{h) is given by 



91 r| 9 /"-f-oo ^—a{m,lc,€) i ( ^ — a{m,k,e)\ri 

e^^{h)=wt''-2\J\ + ^-^ y / dka{m,k,e)-^ , , ' , , , , 

(4.16) 



with 



It is straightforward to check that the /i*(| J|) minimizing e^^{h) for a given value of \ J\ is equal 
to +OC' as soon as \J\ is larger than a critical value | Ji | = J2m I dkae~°'{l + e~°')~^ . This means 
that the minimal energy in-plane state is the staggered antiferromagnet for |J| > |Ji|. If we 
define ef_^^ = e"^^ (h) + 2\J\{1 — h~^), we have ef'^ < e*^, for all /i > 1, and moreover ef^^ < eo, 
with Co defined in (4.13). For h> 1, the critical value of | J| for which e'^-'^(l) = e^^{h) is given 
by 2| J^l = (e*^ — e*^)/(l — h~^). Similarly, the critical value of | J| for which e'^^(l) is equal 
to the specific energy of the out-of-plane staggered state is given by | J°| = (eo — e^^)/2. 

It is now clear that in order to prove (2.8) in the A = +oo case it is enough to show that 
I I > \ Jc\^ for all /i > 1. This simply follows from a computation, whose details we omit. The 
uniqueness property stated in Theorem 3 follows along the same lines as in the ferromagnetic 
case. This concludes the proof of Theorem 3 in the A = +oo case. 



5. THE FERROMAGNETIC CASE: FINITE A AND POSITIVE MASS. 

In this section we prove Theorem 2. We prove it by showing that if e > and A is larger 
then max{Ce — J, 0}, see (2.5), then in the ground state of (2.1) the spins must be either all 
horizontal or all vertical. Note that if this is the case then Theorem 2 simply follows from the 
discussion of previous section. Let us recall that spins now can assume 4 possible directions 
(those parallel to the two in- plane coordinate axis). We shall denote by S_x spin configurations 

in X <Z Wj [in particular we shall denote by spin configurations in A] and by S_ and S_ 
two (arbitrarily chosen) A = +00 infinite volume ground state configurations with horizontal 
and vertical stripes, respectively. 

We need to introduce some definitions. As in the basic Peierls construction we introduce 
the definitions of contours and droplets. Given any configuration S^^, we define Ay = Ay (5^^) 
to be the set of sites at which the spins are vertical. We draw around each x S Ay the 4 
sides of the unit square centered at x and suppress the faces which occur twice: we obtain in 
this way a closed polygon r(Ay) which can be thought as the boundary of Ay. Each face of 
r(Ay) separates a point x G Ay from a point y ^ Ay. Along a vertex of r(Ay) there can be 
either 2 or 4 lines meeting. In the case of 4 lines, we deform slightly the polygon, "chopping 
off" the vertex from the cubes containing a horizontal spin. When this is done r(Ay) splits 
into disconnected polygons 71, ... ,7^ which we shall call contours. Note that, because of the 
choice of periodic boundary conditions, all contours are closed but can possibly wind around 
the torus A. The definition of contours naturally induces a notion of connectedness for the 
spins in Ay : given x, x' € Ay we say that x and x' are connected iff there exists a sequence 
(x ~ xq, xi, . . . , x„ = x') such that x„i, x„i_|_i, m = 0, . . . , n— 1, are nearest neighbors and none 
of the bonds (xm,Xm+i) crosses r(Ay). The maximal connected components 6i of Ay will be 
called V- droplets and the set of ^-droplets of Ay will be denoted by 2?y (Ay) = {5i, . . . ,5s}, 
or simply Vy. Note that the boundaries T{5i) of the V-droplets 5i € Vy are all distinct subsets 
of r(Ay) with the property: ^\^^T{5i) ~ r(Ay). Similarly we can introduce the notion of 
i?-droplets and of the set X'jj(Ajj). A droplet will be either a ^-droplet or an ff-droplet. The 
set of droplets will be denoted by D = U ©y. 
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The same kind of construction allows one to define FM-droplets as the maximal connected 
regions of spins with the same orientation. We shall call FM-contours the boundaries of FM- 
droplets. 

Given the previous definitions, we can now state and prove the main results of this section. 

Lemma 1 (Peierls' estimate - massive case - small J). If X > J + Cs, with 

J, _ fcillogel if £< 1/2 
^ \c2Ee-^ if e> 1/2 

for two suitable constants c\,C2 > (chosen in such a way that in particular > Jo, with 

Jo the constant in (2.6)) then the following is true. Let S_jy be a .spin configuration in A 
and let S €z V be one of its droplets. If S is a V-droplet (resp. H -droplet), the spin config- 

urationJlA coinciding with S_/^ on 6'^ and with S_ (resp. S_ ) on 6 satisfies If (5^) — -^(Xa) > 0- 

Lemma 2 (Peierls' estimate - massive case - large J). Let the same as in Lemma 
1 and A > 0. Then the energy of any spin configuration associated to the set T of FM-contours 
can be bounded below by Efi,[{K) + (J — Cg)'}^^^-^ \^\, where Efm{^) is the energy of the 
ferromagnetic state and |7| is the length of the contour 7. 

Proof of Theorem 2. A consequence of Lemma 1 is that for A > J + the set 
Ay is either empty or the whole A, i.e., there are no contours in the ground state. This, 
together with the discussion in the previous section implies the result of Theorem 2 in the case 
A > J + Ce- Moreover, a consequence of Lemma 2 is that for J > and A > the ground 
state has no FM-contours, i.e., it is ferromagnetic. Since > Jo this in particular means 
that under the same conditions the conclusions of Theorem 2 are valid (see Remark 1 af- 
ter Theorem 2). Theorem 2 follows by the combination of these two results, choosing > 3(7^. 

Proof of Lemma 1 and 2. Let 7 = r(i5) and note that II{S_j^) - H{T_^) is given 
bounded below by (A — J)|7| plus the difference between the self-energies of Sj, and Tj plus 
the difference between the inside-outside interactions of the spins in 5 with the spins in 6'' . By 
the exponential decay of the potential and the fact that Ta coincides on 5 with the A = -l-oo 
ground state, the first difference is bounded below by a positive constant minus, possibly, a 
term of size const. I7I; similarly the second difference is bounded above and below by const. |7|. 
A computation shows that these two constants can be bounded above by c\ log(l/e), for small 
£, and by C2ee~^, for large s, where c\,C2 are two suitable constants. So Lemma 1 is proven. 
The proof of Lemma 2 goes exactly along the same lines: the energy of a state with a non 
trivial set of FM-contours, compared to the energy of the FM state, is given by X]7er 
plus the inside-outside energy associated to any FM-droplet. The latter is bounded below by 
—const. X^-,gr with the constant bounded as discussed above. 



6. THE FERROMAGNETIC CASE: FINITE A AND ZERO MASS. 

In this section we want to discuss how to generalize the Peierls' estimate of Lemma 1 to the 
case e = 0. In this section we shall only consider the case < J < Jo, with Jq defined as in (2.6). 
The main result of this section is a generalization of Lemma 1 to the massless case. Its proof 
requires the use of the screening properties of the A — +00 ground state (related to its striped 
nature, in particular to the fact that its total polarization is vanishing). In order to describe 
the result we also need to introduce the notion of simple droplet: we shall say that a droplet 
5 is simple if either it is simply connected or it winds around the torus and its complement 
is connected. One crucial property of simple droplets 5 we shall need is that the number of 
sites in (5 at a fixed lattice distance d from 5'^ is bounded above by |r((5)|. Note also that any 
collection of droplets T> associated to some spin state S_f^ contains at least one simple droplet. 

We are now ready to state the main result of this section. 

Lemma 3 (Peierls' estimate - massless case). Let £ = 0. IfO<J<Jo and A > const. h*{J), 
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for a suitable constant and with h*{J) the minimizer of the r.h.s. of (2.4), then the following 

is true. Let 5^ be a spin configuration in A and let 6 € V be one of its droplets. If S is a 

simple V-droplet (resp. H-droplet), then the spin configuration T_t^ coinciding with on 6^ 

-.H -> 
and with S_ (resp. S_ ) on 5 satisfies H{S_j^) — H{Tj^) > 0. 

Proof of Theorem 1. An immediate consequence of Lemma 3 is that in the ground state 
I'D] = 1. In fact if by contradiction the ground state configuration had I'D] > 1, then it 
would be possible to reduce the energy by changing S_f^ into in the way described above 
(note that we are using that, as remarked above, any droplet configuration T? always contains 
at least one simple droplet). Then all spins are either horizontal or vertical and this, together 
with the discussion of Section 4, implies the result stated in Theorem 1. 

Proof of Lemma 3. With no loss of generality wo assume that (5 is a ^-droplet. We rewrite 
S_j^ in the form: S_j^ = S_g[J S_gc where the spins in are all vertical and in particular they 
consitute a maximally connected component of vertical spins. We rewrite 

HiSj,) - H{t^) = HoiSg) - Hoisf) + H^iSslSs.) - i?i(^f ||,.) (6.1) 

where Ho{S_g) is the internal energy of the spins in 5 and Hi{§^ \Ssc) is the inside-outside 
interaction between the spins in 6 and those in (5^. 

Let us consider the auxiliary configuration D_/^ — S_gU S_ga coinciding with 5^ inside 6 and 

->v 

coinciding with the (A = +oo) vertical ground state 5^ outside S. (6.1) can be rewritten as 

H{S^) - H{t^) = [H{D^) - H{sl)] + H,{S,\Ss.) + (6.2) 

We shall now separetely estimate the contributions from the three square brackets in the r.h.s. 
By (4.11) and defining 7 = r(^), the first term can be bounded from below by 

H{D^)-H{sl)>-2J\j\+ '«e-"'**/iB+ ^ Ke-"'** (/is, + /isj (6.3) 

where the notation is the same as in (4.11) and we used that cj, dj > Ke~°''^ , see lines following 
(4.11). The second term in (6.2) can be bounded below as 

^^id^l^c) > A|7| - \Hdip{B6\S,.)\- Y \Hd^p{{BuB2}s\Ss.)\ (6.4) 

B:hBjth* <Bi,B2>: 

^B2 

where Hdip{Bs\S_ga) is the dipole-dipole interaction energy between the spins in BDS and 5^c, 
and similarly Hdip{{Bi, B2}s\Rs'=) is the interaction energy between the spins in {Bi U B2} fl 6 
and S_gc . Note that now in the second sum also pairs of blocks with both blocks of optimal size 
are included. The r.h.s. of (6.4) can be bounded below by: 

A|7| - const. [ Y f+ E f]- E \Hdip{{BuB2}5\Ssc)\ (6.5) 

xeB-.hBT^h* ^ xeBiUBa: ^ <Bi,B2>: 

hBi'^^B2—^* hBi=^B2=^* 

where z^^ > 1 is the lattice distance between x and 5'^. In (6.5) we are only left with the 
summation over pairs of optimal blocks. Some of these pairs can be at a distance from 5'^ smaller 
than h*{J), and the contribution from all such pairs is bounded below by —const. h*{J)\^\. Now 
note that a pair of nearest neighbor optimal blocks has vanishing polarization, so we can boiuid 
the contribution of any optimal pair {Bi,B2} at a distance from larger than h*{J) from 
below by 



— const. /i*(J) 'Y^ ~ 



^2 

eBiLIB2: 
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(note that the ^ has been now replaced by 57!)- '^^^ total contribution due to pairs of optimal 
blocks can then be bounded below as: 

-const. /i*(J)|7| - const. /i*(J)^ J2 ^ - -const. /i*(J)|7| 

d>i ^<^>i-- 
- ^^=d 

where the last constant is in general different from those at l.h.s. and we used that the number 
of sites at a distance d from 6'^ is bounded above by I7I (because (5 is a simple droplet). Using 
this result in (6.5) we find that, for a suitable C > 0: 

Hi{Ss\S,.)>{\-Ch*{J))\^\-c[ E f+ E f]- (6-6) 

Similarly, the contribution from the last square bracket in (6.2) is bounded below by 
- C/i* ( J) 1 7 1 . The conclusion is that 

H{S^)-H{Tj,)>{X-W)h\+ E E ('^e-"''*-^) (6.7) 

''Bi?^'*B2=''* 

for suitable constants Ci, C2 > 0. To the purpose of a bound from below we can trow away all 
terms with z^^ > C-^.e"'^ /k in the two sums. We are then left with: 

H{S^)-H{t^)>{\-WM-iC2 E f- (6.8) 

Since 5 is simple, the number of points at a distance d from (5^ is at most |7|. Then the 
last summation is bounded below by —const. |7| log(C2e"'' / k) and the proof of the lemma is 
concluded. 



7. THE ANTIFERROMAGNETIC CASE. 



In this section we conclude the proof of Theorem 3. We already proved it in the A = +00 
case, see Section 4B. So, we now simply need to show that if either (i) \J\ is larger than some 
absolute constant ni and A = or (ii) \ J\ < ki and A larger than some absolute constant K2, 
then the system always prefers to have the spins all oriented either out-of-plane or in-plane 
horizontally or vertically. If this is the case then Theorem 3 follows, with Aq = Ki + K2. 

The easiest case to handle, that wc shall treat first, is the case of \J\ large and A = 0. 
In this case, both the in-plane and the out-of-plane minimal energy states display staggered 
antiferromagnetic order. In analogy with the FM-contours defined in Section 5, we can introduce 
the notion of AF-contours, obtained as union of dual bonds separating nearest neighbor spins 
which are either idential or orthogonal. By construction, the AF-contours separate maximally 
connected regions of spins displaying staggered antiferromagnetic order (AF-droplets) . Let Vi 
(resp. T>o) be the set of in-plane (resp. out-of-plane) AF-droplets and let Cjp = e*^ — 2\J\ 
and Cop = cq — 2\J\ be the in-plane and out-of-plane minimal energies (we recall that e*^ and 
Co were defined in Section 4B). As discussed in Section 4B we have Cip > Cop- Given a spin 
configuration 5^ whose set of droplets is ViUVo, its energy can be bounded below by 



where the two (absolute) constants take into account the dipole interaction energy of the spins 
inside the droplet with the spins outside (note that this inside-outside interaction is simply 
proportional to the length of the contour, because of the screening effect associated with the 
antiferromagnetic phase). By (7.1) we see that if \ J\ is larger than an absolute constant Ki then 
the unique ground state is the out-of-plane staggered antiferromagentic state. 
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Let us now turn to the case |J| < ki and A is larger than some absolute constant K2, to be 
determined below. Given any spin configuration S_j^, let us define ff-droplets and ^-droplets as 
in Section 5 and let us also introduce the notion of 0-droplets as the maximal connected regions 
of out-of-plane spins (the set of O-dropIcts will be denoted by Vq). Since the dipole interaction 
between an in-plane and an out-of-plane spin is zero, the energy H{S_fJ can be rewritten as: 



H{Sa)={^ E \m\+Ho{VHlJVv)) + {^ Yl \m\+HoiVo))^{I) + {n) (7.2) 



where HaiVn U Vy) (resp. H^iVo)) is the interaction energy among the in-plane (resp. out- 
of-plane) spins, corresponding to the Hamiltonian (2.1) with J < and A = 0. Using the 
chessboard estimate, (//) can be easily bomidcd below by {^/'^)^f,e'Do + SdeD„('^'o ^ 

2| J|)|5|. In order to bound (/) we follow the same strategy of Section 6, but we first need to fill 
the regions occupied by the 0-droplets by an auxiliary configuration of in-plane spins. To this 
purpose, we define the auxiliary spin configuration T^, which coincides with S_f^ on Vr U Vy, 

and with S_p^ on Dq ■ We have: 



(7) > 77(f ^) - E [e*\5\ + const.|r(5)|] , (7.3) 



where e* is the specific energy of S_f^ (in terms of the notation introduced in Section 4B, 
e* = e^^{h*) and h* = h*{\J\)). Now H{T_f^) can be bounded below by the method discussed in 
Section 6. However now the role of "blocks of the wrong size" is played by "elementary defects", 
i.e., sequences of 3 contiguous identical spins. Note in fact that none of the minimal energy 
states discussed in Section 4B contains such subconfiguration of spins and a straightforward 
computation shows that the state obtained from an elementary defect by repeated reflections 
has specific energy larger than e^^ {h), for all /i > 1. This makes the present discussion much 
simpler than the discussion of Section 6. The result is that, as long as A is larger than an 
absolute constant, we have that 7f(T^) > |A|e*. Combining this bound with (7.2) and (7.3) we 
finally conclude that 

ff(lA)>(eo-2|J|) E l^l + e* \^\ (7-4) 

and this concludes the proof of (2.8) and of Theorem 3. 



8. CONCLUDING REMARKS. 



We rigorously proved existence of periodic striped order with periods of length h{J) in the 
ground states of a 2D system of dipoles with restricted orientations, with 3D dipole-dipole 
long range interactions competing with a nearest neighbor ferromagnetic exchange interaction 
of strength J. We also considered such system with an antiferromagnetic exchange |J|, in 
which case we proved existence of a reorientation transition from an in-plane to an out-of-plane 
ordered ground state, as \J\ increased. Finally, we gave an example of soft striped order in the 
form of ID sinusoidal spin wave [33]. 

Unfortunately, even the ground states of more realistic models used to describe thin films are 
still far from being solved exactly. On the basis of variational arguments and approximations, 
one expects spontaneous formation of mesoscopic stripes in anisotropic systems of out-of-plane 
spins interacting via the dipole-dipole (or Coulomb) long-range interaction and a short range 
FM exchange, like those considered in [4, 11, 28, 29, 43]. It is also expected that in the presence 
of a sufficiently strong uniform magnetic field, oriented perpendicular to the plane, the ground 
state should exhibit periodic order in the form of bubbles, i.e., domains of spins parallel to 
the external field of quasi circular shape [19]. Striped or bubble patterns arc also expected on 
the basis of an eS'ective (mean field) free energy functional in 2D electron gases [24, 41, 42], 
Langmuir monolayers and liquid crystals [37]. 

Even less is known rigorously for positive temperatures. In particular it is unclear, even 
on a heuristic level, whether the expected striped or bubbled order for discrete spins should 
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have strict long range order (LRO), or rather quasi-long range order (QLRO) characterized by 
order on short scales and a power law decay of the order parameter correlation functions. One 
of the few rigorous results, which we are aware of, about the positive temperature behavior 
of this class of systems is the recent proof by Biskup, Chayes and Kivelson of the absence of 
ferromagnetism in d-dimensional Ising models with long range repulsive interactions, decaying 
as l/r^, d < p < d + 1, interactions [6]. 

The models discussed here are related to a class of systems with Kac potentials {i.e., long 
range potentials of the form 'j'^v{jr)) considered by Lebowitz and Penrose in [27]. There 
they computed the free energy density of a system of particles interacting via a short range 
interaction, favoring phase segregation, and a long range nonnegative definite Kac potential. 
They showed that in the limit 7 — > there is no phase transition in the thermodynamic sense, 
even though the pair distribution function has the form characterizing a phase transition, at 
least over length scales much smaller than 7"^. They concluded that the repulsive Kac potential 
causes the distinct phases of a normal first-order phase transition to break into droplets, or 
froth, of characteristic length large compared to the range of the short range potential and 
small compared to 7"^. We do not know the scale of these domains, or even whether they 
destroy the first-order transition in 2D and 3D when 7 is small but finite. 

In ID we can show, using the method of reference [20], for H = —J'^aiai+i + 
76xp{— 7|? — j\}, that for 7 — > the ground states are periodic with period propor- 
tional to "y~^{J^y^^, i.e., on the macroscopic scale 7"^ the period goes as {JjY''^. This is 
very reminiscent of what happens to the minimizers of continuum energy functionals used 
to model microphase separation of diblock copolymers and many other physical systems 
[2, 10, 12, 13, 32, 34]. 

For Heisenberg spins with long range interactions the Hohenberg-Mermin- Wagner argument 
is not applicable and it could very well be that the long range tails stabilize phases against 
thermal fluctuations, even in two dimensions. This issue has been discussed in some detail for 
the case of rotators in two dimensions, interacting via a 3D pure dipole-dipole interaction [1, 9]: 
this is a case where linear spin- wave theory predicts non-existence of LRO, while non-linear 
corrections and renormalized spin-wave theory seem to suggest that LRO survives at positive 
temperatures. 

For scalar fields describing the local magnetization or electron density in an effective free 
energy functional theory, the general expectation is that in the presence of an anisotropy term 
(possibly induced by the underlying crystalline structure) QLRO should survive at positive 
temperatures. On the contrary even QLRO should generally be destroyed by thermal fluctu- 
ations in the case of isotropic interactions [44]. The issue is rather subtle, however: analyses 
based on a Hartree approximation would generically predict the presence of a first order phase 
transition from a high temperature disordered phase to a low temperature striped phase [7]. 
If both predictions are correct, this class of models would exhibit a rather peculiar first order 
phase transition from a disordered state to a locally ordered state without strict LRO. The issue 
has been investigated by Jamei, Kivelson and Spivak [24, 41, 42] in the context of 2D electron 
gases and by Tarjus, Kivelson, Nussinov and Viot [44] in the context of the frustration-based 
appoach to the glass transition in supercooled liquids and structural glasses. Jamei, Kivelson 
and Spivak argued that for 2D electron gases, in the presence of Coulomb interactions, first 
order transitions (as those predicted by Brazovskii) are not allowed and a sequence of transi- 
tion between different mesoscopic patterned states should generically appear instead. On the 
contrary, as discussed in [44], numerical simulations of frustrated spin models with Coulomb 
interaction predict a finite jump in the energy at a critical line separating the paramagnetic 
state from a locally order striped state. 

Refinements of the mean field or variational arguments seem very difficult: the natural ef- 
fective continuum theories describing systems whose interaction has a Fourier transform with 
a non-trivial minimum in fc-space are, at least naively, non-renormalizable [23, 40]. It would 
be very interesting to provide convincing arguments for the patterned states to be (globally) 
stable against the presence of thermal or quantum fluctuations, as well as against "generic" 
perturbations in the form of long or short range interactions. 
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APPENDIX A: A NEW PROOF OF PERIODIC ORDER IN 
REFLECTION-POSITIVE ID ISING SYSTEMS 

We mentioned in the Introduction that the proofs of Theorems 1 and 2 in [20] contain a 

technical error, which was caused by overlooking some exponentially small terms that came 
from the choice of periodic boundary conditions. In this Appendix we show how to repair this 
error. We exploit a new method, which is simpler than the one used in Section III of [20] 
and which docs not make use of periodic, but rather of open boundary conditions, as will be 
discussed below. The proof given in this Appendix also implies the results discussed in Section 
4. 

First of all, let us point out the mistake in [20] : a term 

(-l)i-J(l-e-«''')(l-e-"'^Oe"^"^ H (^"^^ 

l<i<j<M/2 i<k<j 

is missing in the definition of Hfj{a,hj^) in (1.9), i.e., in Eq.(9) of [20]. Taking into account 
this term, one can check that, whenever the total length Nl = ^M^i<ohi of the blocks 
on the left is different from the total length A^^ = J2i<i<M. hi of the blocks on the right, 

the new energy functions H{a,{h]^,6hi)) ^'^^ H{a,{6hji,hj^) obtained after the first re- 
flection, see Eq.(ll) of [20], are not periodic. For this reason it does not seem possible to 
repeatedly reflect, as one must do in order to obtain the checkerboard estimate, see [20], Eq.(12). 

Let us now show how to correct the proofs of Theorems 1 and 2 in [20]. Let us denote Eq.(x) 
of [20] by (I.x) and let Hm be the periodic Hamiltonian in (1.1): 

N 1 
HN{a) = - J ^ (Jiai+x + - ^) '^i ' "^f - ^) = 5Z + 2nN\P 

i=-N+l -N+l<i<j<N „ez<* I J "T I 

(A.2) 

where aN+i = ct-n+i and p > 1. As discussed in [20], Theorems 1 and 2 are consequences of 
the chessboard estimate, used in [20] and stating that 

n n 

H{h,,...,hn)>Y,hie{hi) , Yhi = N, (A.3) 

i=l i=l 

where H{hi, . . . ,hn) is the energy H]sr{a) of a spin configuration whose corresponding block 
configuration is h = {hi, . . . , and e(/i) is the energy per site of the infinite system with a 
periodic configuration of blocks all of the same size and alternating sign (let us recall that a 
block is a maximal sequence of spins all of the same sign) . In the following we shall also need to 
introduce the analogue of H{h) with open boundary conditions, to be denoted by H°{h). Since 
the long range potential is summable, we have 

mtimes 

H{h) = Hm — F°(£'^"?^) (A.4) 

m^oo m 



and we find that the chessboard estimate (A.3) is a consequence of the following: 
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Chessboard estimate with open boundary conditions. Given a finite sequence of 
blocks -Aj = {h\,. . . , hl^.}, let us denote by 9A^ = {hl^., . . . ,h\} the reflection of (with the 
sign of the spins in the reflected blocks being opposite to what they were originally). By e{A^) 
we denote the infinite volume energy per site of the configuroMon (. . . , A^, OA^^Aj^, OA^, . . .) and 

let ttj = Y^Y^i hj < oo. Then, for any collection Aq, . . . ,A^^i, with n > 1 and Xir=5)^ '^i ~ ^> 
we have 

n 

H°{Ao,Ai,...,A„,A„^.^) > (ao + a„+i)eo + ^a„e(A„) , (A.5) 

where cq is the infinite volume specific ground state energy. 

Proof. Let us first note that, for any sequence of blocks A, e{A) > aeo, where a = J2heA ^• 
This can be proven as follows. Denote by the spin configuration corresponding to A, so that 
H'^iA) = Hl{n^). If ^ = ±, we have H^,{a^,^a^) - 2ffO(aJ+^S,„t(^„;aJ, where K,^'aJ 
is the spin configuration of length 2a obtained by juxtaposing and uia^, and ivEintis^', 3^) is 
the interaction energy between the two halves. For one of the two choices ui = ± this interaction 
energy is nonpositive, so that 2H^{a^) > mmg_H2a{a). Iterating, we find: 

H^ia,) > hm 2-™ mm H^^,{a) = aeo , (A.6) 

m—>(x> a_ 

which is the desired estimate. 

Next, let us recall that reflection positivity of the long range potential l/r^ implies the 
following basic estimate. 

Lemma A.l. Given two finite sequences of blocks A_ = {h-M+i, ■ ■ ■ ,ho} and A_^ = 
{hi, . . . , hN}, with M, N > 1, let 0A_ = {ho, . . . , /i_m+i} and 6A_^ = {Hn, ■ ■ ■ , hi} be their 
reflections. Then we have: 

H°{A_,A^) > ^H%eA^,A^) + ^H°{A_,eA_) . (A.7) 

Now we are ready to prove (A.5). We proceed by induction. If n = 1 in (A.5), then by 
reflection positivity, i.e., by Lemma A.l, we have: 

H%Ao,A„A2) > A2) + lH%Ao,A„eA„eAo) (A.8) 

By (A.6), the first term in the r.h.s. can be bounded from below by 0260- The second term can 
be bounded by a second reflection: 

lH°{A„A„0A„eA,) > i7J°(Mo,A) + Jff°(^o> Ui)^', ^^0) (A.9) 

where by definition (AJ®'' = (A^, M^, A^, M^). By (A.6), the first term in the r.h.s. of (A.9) 
can be bounded from below by aoeo/2, so we end up with: 

if°(^,4i,42) > «2eo + iaoeo + Ji?°(Ao, {A,f\eAo) (A.IO) 

Iterating we find: 

H°{Ao,A^,A2)>a2eo + aoeoy2-''+ lim 2-"H°{Ao,{A^f^" ,0Ao) (A.ll) 

n>l 

Note that the last term is equal to aie(Aj), so the desired bound is proven for n = 1. 

Let us now assimic by induction that the bound is valid for all 1 < fc < n — 1, n > 2, and let 
us prove it for fc = n. There are two cases. 



19 



(a) n = 2p for some p > 1. If we reflect once we get: 

iAp+1, • • • )di2p+l 

If we now regard A!p+i = i^Ap+iiAp+i) and -A^ = {A^, 0Ap) as two new sequences of blocks, 
the two terms in the r.h.s. of (A.12) can be regarded as two terms with n = 2p — 1 and, by the 
induction assumption, they satisfy the bounds: 

2p 

^^°(^^p+i, • • • • • -'42^+1) > 2a2p+ieo + 2 ^ a^e^.) (A. 13) 

i=p+l 

P 

^?°(^,^i, • • • ,4' • • -i^AiiOAo) > 2aoeo + 2^aie(A.) (A. 14) 

1=1 

where we used that, by construction, = 2ap, ap_^_i = 2ap+i, e(-Ap) = e{Ap) and e{A^p_^_i) = 
e(-Ap_|_i). Therefore, the desired bound is proven. 

(b) n = 2p+l for some p>l. If we reflect once we get: 

H'^iAoiAu ■ ■ ■ ):d2p+l)ii2p+2) — (A. 15) 

-if'^(M2p_|_2) • • • ) 0Ap+2tAp+2i ■ ■ ■ 1A2P+2) + 2-'^'^(— 0'— 1' • • • )4+l) ^4+1) • • • ) ^Al,^Ao) 

The first term in the r.h.s. corresponds to n = 2p so by the induction hypothesis it is bounded 
below by a2p+2eo + X^i=p+2 ^i^(— i)- As regards the second term, using reflection positivity 
again, we can bound it from below by 

\h\A„,A„ ...,Ap,0Ap,..., Ml, Mo) + ^ff°(Ao, Ai, ...,Ap, {Ap+i)^\ OA^, . . . , Mi, Mq) 

(A.16) 

By the induction hypothesis, the first term is bounded below by aoeo/2 + (1/2) J2^=i ^^^(-AJ, 
and the second can be bounded using reflection positivity again. Iterating we find: 



-^^''(ilo'^l' • • • '^2p+l)^2p+2) ^ (A-17) 
2p+l p 

> a2p+2eo + Yl (^i^iAi) +{Y1 2"") ("0^0 + £ "'i^iAi)) + 

i—p-\-2 n>l i—1 

hm 2-"i?oUo, Ai, ...,Ap, (Ap+i)«2", Mp, . . . , Mi, 9^,) . 



n— »oo 



Note that the last term is equal to ap+ie(Ap^i), so (A. 17) is the desired bound. This concludes 
the proof of the chessboard estimate with open boundary conditions and, as mentioned above, 
of (A.3) and of Theorems 1 and 2 in [20]. ■ 
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